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Abstract 



We present upgraded theoretical predictions for inclusive and ex- 
clusive radiative decays of B mesons. Our results include those next- 
to-leading order corrections that have already been computed. Our 
^ ' best estimates in the Standard Model are BR{B K*'j) = (4.3 it 

^ '. 0.9l{:^) X 10-^ BR{B Xs-f) = (1.9 ± 0.2 ± 0.5) x lO^^^ T{B 

K*i)/T{B Xs-f) = 0.23 ± 0.09 and BR{B X^g) = (1.57 ± 
O.lS+olg ± 0.23) X 10-3. We also consider limits found with two- 
Higgs-doublet models. 
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Introduction 



Radiative decays of B mesons represent very important tests of the weak 
interactions and of the role of effective flavour-changing neutral currents. 
Among these decays, 6 — 57 and b ^ sg are theoretically clean and sensitive 
to physics beyond the Standard Model, e.g. charged scalar Higgs models 
and/or SUSY models The experimental measurement of the exclusive 
branching fraction BR{B K*-f) = (4.5 ± 1.5 ± 0.9) x 10"^ by the CLEOII 
collaboration [^] and the imminent measurement of the inclusive rate, on 
which an upper limit BR{B -> X,7) < 5.4 x 10"^ (95% C.L.) [| already 
exists, offer the opportunity to compare experimental results and theoretical 
predictions for these quantities. 

From the theoretical point of view, the prediction of the rates consists as 
usual of two steps. 

On the one hand, it is necessary to compute the renormalization of the 
coefficients of the effective Hamiltonian to take into account the effects of 
strong interactions at short distances, i.e. for scales rrih < ft < My/. It turns 
out that renormalization effects at the leading order (LO) have important 
consequences, since they almost double the amplitude obtained without their 
inclusion ^. Unfortunately, the full set of next-to-leading corrections 
to these decays, which are necessary for a consistent use of Aqcd, of the 
renormalization scale, and for more accurate predictions, are not available 
yet. 

On the other hand, it is necessary to compute the hadronic matrix ele- 
ments of the operators appearing in the effective Hamiltonian. For inclusive 
decays, in the framework of the Heavy Quark Effective Theory (HQET), 
it is possible to predict the rate using the parton model, with computable 
corrections that are expected to be of order I /ml 0-|n|. For exclusive 
decays, one has to know the relevant hadronic form factors, obtained from 
a non-perturbative estimate. For B K*'~f there is only one form factor, 
which we will denote in the following by Fi{0). The exclusive channel had a 
bad reputation because different predictions of the rate varied by orders of 
magnitude. In the recent past, however, lattice QCD [|l^, [l^ and QCD sum 
rules |T^-[0 have procured more reliable results and the theoretical uncer- 



tainties have been substantially reduced. This makes the exclusive decays 
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more interesting as tests of the Standard Model. 

We present an upgraded analysis of the inclusive and exclusive 6 — *■ S7 
and of the inclusive b ^ sg decay rates, which takes into account several 
improvements made recently: 



The coefficient of the magnetic (chromo-magnetic) operator, C7 (Cg) 
is now established. In refs. [0, |18[, the long-standing problem of 
the regularization dependence of the LO coefficients of these operators 
was solved; the results were found to be different from all previous 



calculations |jT9[. Those of refs. [0, |18[ were subsequently confirmed 
inref. 



The next-to-leading order corrections to the anomalous dimension ma- 



trix are partially known |2J] . We will make use of this information 
to try to evaluate the effect of the next-to-leading terms and include it 
in the error on the final predictions. We also compare results obtained 
in the HV and NDR regularization schemes, since a spurious regular- 
ization dependence is introduced by the incomplete NLO terms. The 
NLO corrections that we can already include diminish C7 in both the 
HV and NDR schemes. It turns out that in HV the dependence on 
the renormalization scale is substantially reduced. A similar effect is 
also found in NDR, where nevertheless a sizeable dependence is still 
present. We will also make use of the partial calculation of the 0{as) 
corrections to the inclusive rate, which has been computed in ref. [p^] . 

We will study both the inclusive and exclusive channels, the latter being 
ignored in most of more recent analyses. We will make use of lattice 
and QCD sum rules predictions for the relevant form factor -Fi(O). 

By varying within their errors the experimental and theoretical quan- 
tities, we obtain a distribution of values for the theoretical predictions, 
from which we estimate the theoretical uncertainty. 



Our predictions contain several differences with respect to the recent 
analysis of ref. [^]: 

i) the estimates we present are also for the exclusive B — > K*'y rate; 
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ii) according to refs. |29[, we have used the running mass, and not 
the pole mass, for the evaluation of the inclusive rate; 

iii) 0(1/ ml) terms have been taken into account in the present study; 



iv) NLO corrections to the coefficient function ^Tj and to the anomalous 
dimension matrix T2\ have been included by us in evaluation of the values 
and theoretical uncertainties for the rates. 

Of all these differences, the most important is the last one, since the 
known NLO corrections diminish the strong fi dependence and reduce the 
values of the rates. Apart from this, the bulk of our results is substantially 



in agreement with those of ref. p3 . 



Basic Formulae 



The effective Hamiltonian responsible for 6 — 57 and b sg decays can be 
written as §, §, 

where: 

Qde 



O7 
0« 



167r2 "V )r 



and is the renormalization scale of the relevant operators. From the ef- 
fective Hamiltonian one can derive the decay rates. In the following, we 
report the main formulae that will be used in the numerical evaluation of the 
inclusive and exclusive branching fractions. 

Inclusive B Xg'-^: 

-I th 



BR{B ^ X,7) 



T{B ^ X,7) 
T{B ^ Xlui) 



X BW''{B Xlui), 



where : 



r(g-^x,7) 

T{B Xlui) 



-[ th 



\Vcb? Qirgimjmb) 



X Fx \C^^^{fi)\ 



(2) 



(3) 
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with the phase-space factor g{z) given by: 



24z^ \n{z) 



(4) 



and 



K{mt/Mw,l-i) 



(5) 



fi(mc/mb,/i) 

In eq. (|^) the quantity fl{z) contains the 0{as) QCD corrections to the 
semileptonic decay rate Within a good approximation it is given by 

2a^(/i) 



n{z, /i) ^ 1 



TT 



3 

+ 2 



(6) 



the factor K{mt/M^, fi) contains the 0{as) NLO corrections to the B Xg'j 
rate, due to real and virtual gluon emission ||2^. The calculation of ref. [27 
does not contain the full set of 0{as) corrections, which would require a two- 
loop calculation of the coefficient function and a three loop calculation of the 
anomalous dimension. In ref. the authors only included some terms, 

which one may argue to be the most important ones. We have used the 
results of ref. |27| in our analysis, in the same spirit in which we considered 



the NLO anomalous dimension, which is only partially known. 

The scale n, which appears in some of the previous formulae, is there for 
two reasons. In the numerator of eq. (|^), /i denotes the renormalization scale 
of the effective 6 — *• 57 Hamiltonian. It also represents the scale at which we 
decided to compute the expansion parameter as for the QCD corrections to 
the semileptonic decay rate. For simplicity we have taken the same value of 
/i in both cases. According to refs. p8|, in the framework of the HQET, 
to avoid problems with renormalon (non-perturbative) effects, we have used 
as expansion parameter the running mass rrih = mb{n). For this reason, since 
the b ^ s'y operators are evaluated at a generic /i, we have omitted the factor 
{mb{fi) / mb{pole))'^ which was included in F by the authors of ref. ||23l 
eq. (D. 

Exclusive B K*'j: 



see 



BR{B K*-f) 



-1 th 



T{B ^ K*-f) 
TiB ^ X,7) 
BR^'^P^B Xlvi) 



X 



T{B ^ X,7) 



-[ th 



T{B ^ Xlui) 



(7) 
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where 



th 



rrib 



1 - 



Ml, 



l + (Ai-9A2)/(2m2; 



xlFi 



(8) 

In the HQET formahsm, the parameters Ai and A2 describe the leading non- 
perturbative corrections (at order 0{l/ml)) to the parton model predictions 
for the inclusive rate [|ll|]. They are related to the kinetic energy of the 
h quark (inside the B meson) and to the B-B* mass splitting. The 
corrections cancel in r{B —>■ Xs'y)/r{B —>■ Xlvi) but not in the ratio (|^); 
Fi(0) is the relevant form factor defined by: 



{K\p\r,)\sa,,{l + ^,)q^h\B{p)) 



+ 2 



V:iMl-Ml,)-iv*-q)ip + p'), G,{q 



(9) 
(10) 



with ^2(0) = Fi(0)/2. 
Inclusive B XsQ: 

The theoretical estimate of BR{B X^g) is at present much rougher 
than BR{B —>■ Xsj). Even a partial evaluation of the coefficient function and 
1/ml terms is missing in this case. We give our prediction, which improves 
with respect to previous calculations only because Cg'^^ has been corrected 
and we include the NLO terms of the (6 x 6) anomalous-dimension submatrix. 
To evaluate the branching ratio we have used the formula: 



BR{B Xs 



r{B Xs 



T{B ^ Xlui) 



th 



X BR^''P(B Xlui), 



fill 



where : 
T{B 



TiB^Xlui)^ 



th 



\vzvtt\' 



2as(yu) 



x|C3^^/(^)p. (12) 
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Parameters 


Values 




0.95 ±0.04 


rric/mb 


0.316 ±0.013 


mt(GeV) 


174 ± 16 


Ai(GeV") 


-0.15 ±0.15 


A2(GeV2) 


0.12 ±0.01 


mb(/i = m;,)(GeV) 


4.65 ±0.15 


Fi(0) 


0.35 ±0.05 


BR^'-'P{B Xlvi) 


0.107 ±0.005 


AS^CD(GeV) 


0.240 ± 0.090 


^^ 


mfc/2-2mb 



Table 1: Values of the different parameters used to predict the inclusive and 
exclusive radiative b decay rates. 



Predictions and Uncertainties 
for the Decay Rates 

We give in table |l| the range of variation of all the quantities appearing in 
eqs. which have been used to obtain our results. 

Using the central values of table jT], we show in fig. |l| the fi dependence 
of Cj^^{n) at the LO and NLO. For LO we mean that we have taken the 
anomalous dimension matrix at the LO. For the total rate, we put to zero 
all the 0{as) terms that appear in eqs. including those relative to 

the semileptonic rate. In this case one has also to use K{mt/Mw,f^) = 1- 
In the NLO case we turn on all the known NLO corrections, including the 



(6 X 6) anomalous-dimension matrix computed in refs. ||2T|, In this 

case we have varied K{mt/Mw, fJ') according to the results of ref. i.e. 
0.79 < K{mt/Mw,lj) < 0.86 for mb/2 < /i < 2mb. We notice that the fi 
dependence is reduced both in the HV and in the NDR cases. If we call 
R = {C'r"i^^)\,=m,/2/C'7"i^^)\,=2m,)^ we get ~ 1.72, R^^^ ~ 1.54 and 
R^^ ~ 1.25 []. One would prefer HV because of the reduced fi dependence, 

^This means that the coefficient itself varies only of ^ 10 — 15% in a range of scales as 
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BR{B Xs-f) X 10^ 


(GeV) 


LO 


NLO//y 


NLOjvDi? 


mb/2 


3.81 ±0.47 


1.92 ±0.19 


2.77 ±0.32 


nib 


2.93 ±0.33 


1.71 ±0.18 


2.25 ±0.25 


2nib 


2.30 ±0.26 


1.56 ±0.17 


1.91 ±0.21 




BR{B K*-i) X 10^ 


;u(GeV) 


LO 




NLOjvDi? 


nib/ 2 


6.9 ±1.5 


4.4 ±0.8 


6.4 ± 1.3 


nib 


5.3 ±1.1 


3.8 ±0.8 


5.0 ± 1.0 


2mb 


4.2 ±0.9 


3.3 ±0.7 


4.1 ±0.8 




BR{B Xsg) X 10^ 


//(GeV) 


LO 




NLOatdr 


mb/2 


4.21 ± 0.38 ± 1.03 


1.81 ±0.17±0.14 


3.25 ±0.30 ±0.66 


nib 


2.66 ±0.24 ±0.47 


1.45 ±0.14 ±0.10 


2.19 ±0.21 ±0.32 


2mb 


1.81 ±0.17± 0.25 


1.19 ±0.12 ±0.09 


1.58 ±0.15 ±0.19 



Table 2: Theoretical predictions for exclusive and inclusive radiative B de- 
cays. The last error on BR{B Xsq) comes from Aqcd- We have preferred 
to show it separately, since the rate is directly proportional to a^. 



as shown fig. |I|. In the absence of a complete calculation, however, we cannot 
decide which of the results, NDR or HV, is closer to the real NLO result. 
For this reason, for all the quantities reported below, we will combine the 
results obtained with NDR and HV with mb/2 < fi < 2mb, and include the 
differences in the final estimate of the error. As central value we will take the 
average between the NDR and HV result. Regarding all the other quantities, 
we allow them to vary within the ranges reported in table |^, with Gaussian 
distributions for the experimental parameters and fiat distributions for the 
theoretical ones. In a given regularization scheme and for a fixed value of /i, 
this procedure generates a pseudo-Gaussian distribution of values (see fig. |^) 
from which we deduce the error. For more details, see ref . . 

We want to add some comments on the value and error of -Fi(O) which 

large as ^ ~ 2 — 9 GeV. 
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has been used in the exclusive case. Recent lattice and QCD sum rules 
calculations of this quantity give: Fi(0) = 0.20 ± 0.02 ± 0.06 p|, Fi(0) = 
m and Fi(0) 



o.3o;^° 



0.35 ± 0.05 0, Fi(0) = 0.32 ± 0.05 [|T§, Fi(0) 



16 . For this reason we have chosen to use 



0.310 ± 0.013 ± 0.033 ± 0.060 
-Fi(O) = 0.35 ± 0.05, which covers most of the theoretical predictions (see 
table p. From the numbers given in table |] we quote BR{B —>■ K*^) = 
(4.3 ± 0.91};^) X 10-^ BR{B X,7) = (1.9 ± 0.2 ± 0.5) x 10"^ and T{B 
K*'~f)/r{B Xs'~f) = 0.23±0.09. The first error comes from the width of the 
pseudo-Gaussian distribution of the theoretical values, the second includes 
the /i dependence and regularization dependence of the results. 

The results for BR{B X^g) are reported in table ^ Because of the 
factor as(/i) present in eq. ([T^), the dependence on Aqcd is stronger than 
for BR{B —>■ Xsj), and we have reported the corresponding uncertainty 
separately in the table. Our best estimate is BR{B Xsq) = (1.57 ± 
0.15to;59 ± 0.23) x 10"^ where the last error is due to Aqcd- 



Constraints on 2H models 
from radiative B decays 

We consider the popular 2H model known in the literature as Model II [|^. 
We look for bounds in the plane Mh+ — tan /3, by imposing the experimental 
constraints on the branching ratios for B — > Xg'y and B K*'y. We modify 
the initial conditions of Cj^^ {My/) to take into account the charged Higgs 
contributions ||31[] and calculate eqs. (0)^(0) and (0)^(11)) by varying the 



parameters with the same criteria as before. In presence of charged Higgs, 
the corrections of 0{as) to Ci{Mw), i = 1, ... ,6 are known and were 
included at the NLO in the present calculation. We look for bounds by 
varying the parameters MH+,ta.nj3 in the range Mh+ > 91 GeV, tan/5 < 
0.48 GeV-^xMH+ O-H- The first condition comes from the lower bounc. 



on the Higgs mass established from the limits on Higgs production at LEP. 
The second is obtained from the comparison between theoretical calculations 



and the experimental measurements of BR^^'^^B — > Xrz/^) [Q. Since 
the exclusion regions are substantially independent of tan [3 for large values 
of this parameter, the results will only be given for tan/3 < 3. We observe 



8 



that the exclusion region is strongly reduced by going from the LO to the 
NLO. This is mainly due to the sizeable reduction of the rate when the known 
NLO corrections are considered. Because of the larger theoretical error, the 
exclusive channel cannot be used to put a limit. The results are given in 
fig. H. In view of future measurements of BR{B — > Xs'j), and biased by our 
theoretical predictions, we also give the exclusion region corresponding to an 
upper limit for BR{B X,7) of 4 x 10-^ 3 x 10-^ and 2 x 10"^. Finally, 
in fig. ^, BR{B Xg'y) is reported as a function of Mh+ for tan/3 = 2. In 
the figure, the band indicates the theoretical uncertainty. 
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Figure 1: C7^^(/x) as a function of fi in the LO, NLO-HV and NLO-NDR 
cases. For the definition of NLO see the teoct. 



13 



600 - 



500 



400 



300 - 



Entries 9000 
Mean 0.1707E-03 
RMS 0,1771E-04 



200 - 



100 







0.1 0.12 0.14 0.16 0.15 

BR(B^X,7) 




0.2 0.22 0.24 

X 10" 



Figure 2: Distribution of values for BR{B — > Xg^y) for /i — m^. 
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tan /? - M 
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Figure 3: Exclusion regions (on the left of the figures) in the plane Mh+ — 
tan/3, at next-to-leading order, for different upper bounds of BR{B — > Xsj). 
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Figure 4: Predictions for BR{B ^^7) in the 2H model with tan (3 — 2 are 
given as a function of Mh+ ■ 
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